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PREFACE 


RAND,  among  other  organizations,  pioneered  in  the  use  of  the 
Monte  Carlo  method  for  the  solution  of  con^ilicated  problems  involving 
the  transport  of  nuclear  radiation  such  as  gamma  rays  and  neutrons. 
This  method  has  been  useful  chiefly  when  solution  by  other  means 
was  prohibitively  conqilicated. 

One  of  the  main  disadvantages  of  the  Monte  Carlo  method  or, 
for  that  matter,  any  other  method  of  random  sampling  is  that  one 
never  quite  knows  the  accuracy  of  the  answers.  In  this  report  we 
solve  by  means  of  Monte  Carlo  a  problem  with  a  known  analytical 
solution.  This  problem  is  sijnilar  to  those  occurring  in  actual  RAND 
Monte  Carlo  codes.  Comparison  of  the  Monte  Carlo  solutions  with  the 
known  solutions  allows  us  to  get  a  better  feeling  for  the  accuracy 
of  Monte  Carlo  results  in  practical  problems. 
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SUMMARY 


It  is  shown  that  the  probability  density  function 


dp 


2dx 

x3 


X  >  1 


is  a  reasonable  approximation  to  sampling  distributions  found  in 
several  Monte  Carlo  transport  codes, 

A  machine  code  was  set  up  which  was  used  to  obtain  10,000  runs 
of  10  samples  each  from  the  above  distribution  function.  Histograms 
of  the  probability  density  function  for  the  mean  value,  of  the 
standard  deviation  of  the  mean  and  of  Student*s  ratio  of  10  samples 
were  obtained.  Because  the  standard  deviation  has  an  infinite 
expected  value,  it  is  of  doubtful  usefulness  in  measuring  the 
probable  error  of,  or  placing  confidence  limits  on,  the  mean  value 
of  the  sample. 

The  results  obtained  in  this  paper,  however,  show  that  with 
proper  precautions  the  sarple  standard  deviation  may  be  used  to 
establish  confidence  limits  on  the  sample  mean  in  spite  of  the  fact 
that  it  has  an  infinite  expected  value. 

A  further  examination  of  the  results  of  this  problem  sheds  light 
on  other  properties  of  Monte  Carlo  solutions,  including  the  fact 
that  the  results  are  biased.  Although  the  results  of  this  problem  do 
not  allow  one  to  jsake  quantitative  improvements  in  the  interpretation 
of  Monte  Cjirlo  results,  they  do  uncover  several  qualitative  features 
of  the  Monte  Carlo  results  that  allow  for  more  lucid  interpretation. 
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I.  imODUCTION 

Over  the  past  several  years  questions  have  continuously  arisen 
as  to  the  accuracy  of  solutions  obtained  fron  various  transport  codes 
that  use  the  Monte  Carlo  nethod  of  solution.  The  first  attempts  to 
delineate  the  accuracy  made  use  of  the  calciilation  of  the  sairple 
variance  of  the  sample  mean  as  weU  as  the  sample  mean  itself  for 
each  point  in  the  sample  space.  All  present  RAIID  Konte  Carlo 
transport  codes  calculate  the  sample  variance  of  the  sample  mean, 
and  this  quantity  is  used  to  estimate  the  error  in  the  sample  mean. 
However,  it  has  been  realized  that  this  estimate  of  the  error  may 
oftentimes  be  misleading. If  the  distribution  function  being 
sampled  is  Gaussian  or  a  reasonable  facsimile  thereof,  then  the 
sample  variance,  or  more  correctly  the  sample  standard  deviation, 
iidll  indeed  give  a  Icnov/n  measure  of  the  error.  If  the  distribution 
being  sampled  deviates  from  a  Gaussian  distribution,  the  estimate 
of  the  error  -which  can  be  inferred  from  the  sample  standard  deviation 
of  the  sa'.plo  mean  becomes  more  muddled.  In  particular,  vdien  the 
distribution  function  being  sampled  i:as  long  tails,  the  difficulty 
mentioned  above  ma;/  become  quite  aggravated. 

In  tias  report  wc  will  set  up  a  distribution  f  notion  -wl-dch  is 
sii-rilar  to  ones  enco'untered  in  the  RiM-ID  neutron  and  gam-na  ray 
transport  codes,  i'o'wever,  -with  simplifications,  this  distribution 
function  will  have  a  known  or  calculable  expected  value.  The  Monte 
Carlo  method  is  then  used  to  obtain  estimates  of  this  expected 
value.  The  sample  standard  deviations  are  also  estimated,  and  their 
value  in  estlmting  th.e  error  is  then  exa’raned. 


II.  APPROXIMATION  TO  THE  SOLUTIONS  OBTAINED  HY  >SDNrE  CART.n 


The  RAND  neutron  and  garnia  ray  transport  codes  have  been  used 
to  obtain  Monte  Carlo  solutions  for  a  large  nvuriber  of  problems  wliich 
involve  point  sources  of  neutrons  or  gamma  rays.  All  of  the  solutions 
obtained  ivith  tnese  ^tonte  Carlo  codes  have,  in  general,  tiie  properties 
that  at  close  distances  to  the  source,  the  solutions  have  small 
fluctuations  from  one  range  interval  to  tlie  ne;:t.  These 
fluctuations  increase  as  the  range  increases,  finally  becoming  so 
very  large  that  the  solutions  are  usefiiL  only  out  to  a  certain 
maxiraum  range,  v/iuch  is  commonly  of  the  order  of  10  moan  free  paths 
of  the  radiation  being  considered.  Oftentimes  results  are  needed  or 
desired  at  the  very  rmcLmum  ranges  at  v/iiich  the  Ilonte  Carlo  solutions 
are  usefij.,  and  various  so-called  tricks  are  used  to  make  deeper 
penetrations  possible. 

Even  by  these  methods  it  commonly  happens  tliat  the  fluctuations 
statistical  uncertainties  in  the  solutions  at  large  ranges  arc 
sizeable.  Since  in  many  cases  it  is  necessary  to  live  with  these 
large  errors,  one  woiild  like  to  interpret  then  as  precisely'  as 
possible. 

For  any  particular  element  of  the  sampling  space,  the  Monte  Carlo 
estimates  a  sai:ple  mean  in  accordance  v/ith  the  foUo^ving  formula: 

_  Ex. 

X  -  "if  »  (1) 

vfhere  the  x^  are  the  individual  saTple  values  of  the  quantity  being 
estimated  and  II  is  the  total  size  of  the  sample.  The  sample  variance 
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of  the  san^xLe  mean  is  estimated  at  the  same  time  by  means  of 


The  sanple  standard  deviation  is  simply  the  square  root  of  this 
quantity.  Ordinarily,  if  one  assumes  that  the  distribution  function 
of  is  Gaussian,  then  the  probability  that  x,  the  sa'iple  estimate 
of  the  mean,  vdll  differ  from  the  true  mean  by  more  than  one  standard 
deviation  in  either  direction  is  about  0,32,  and  that  it  vrill  differ 


by  more  than  various  other  proportional  parts  of  a  standard  deviation 
is  in  accordance  with  the  cumnulative  value  of  a  standard  Gaussian 
(normal  0,l).^ 

V/e  vdll  nov;  proceed  to  examine  the  validity  of  such  assumptions 
for  a  particular  distribution  function,  namely. 


X  >  1 

>  1. 


(3) 


This  distribution  function  is  a  reasonably  appro;d-mate  representation 
of  the  distribution  functions  encountered  in  the  R/U®  Monte  Carlo 
transport  codes  as  shovm  in  the  Appendic:.  The  e;qDected  or  mean  value 
of  X  in  Sq.  (3)  is  obviously  2.0.  On  the  other  hand,  the  second 
mom.ent  of  x,  and  hence  the  variance,  is  infinite.  Since  the  variance 
is  infinite,  this  imediately  tlirows  doubt  on  the  whole  idea  of  usin^ 
the  sample  variance  as  com.puted  ty  means  of  Eq.  (2),  since,  if  enough 
samples  are  talcen,  it  must  become  infinite. 

^:-The  use  of  the  standard  deviation  of  the  sample  mean  for 
establishing  confidence  limits  on  the  sample  mean  is  discussed  again 
on  p.  11  and  in  Ref.  .,  pp.  517-519. 
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III.  MOmE  CARLO  SOLUTION  OF  PROBLEM 


A  machine  code  ;vas  set  up  to  obtain  the  solution  by  reindon 
sampling  for  the  mean  value  of  x  as  given  by  the  distribution  function 
of  Eq.  (3) .  Samples  from  the  distribution  function  vrere  obtained  in 
one  of  the  usual  ways  setting  the  cummulative  distribution  function 
equal  to  a  random  number  uniXoinly  distributed  between  0  and  1. 
Symbolicallj’- 


P  = 


R, 


(4) 


or 


X 


1_ 

yR* 


One  then  obtains  x  for  a  sample  size  N  by 


(5) 


X 


1 

N 


(6) 


The  sai^ple  variance  of  the  sample  mean  is  given  by 


— o 

o'' 


ir’-  i=i  \ 


(7) 


It  should  be  noted  that  the  variance  given  in  the  above  equation  is 
not  the  sample  variance  itself,  v/hich  v/ould  in  fact  be  II  times  the 
above  quantity,  it  is  rather  the  sample  variance  of  the  sample  mnan; 
hence  presumably'  a  measure  of  the  statistical  error  in  the  quantity 
X.  Besides  computing  the  sample  standard  deviation  of  x,  vrfiich  is 
simply  the  square  root  of  the  quantity  given  by  Eq.  (7),  the  code 
also  computed  Student’s  ratio 
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The  nunber  of  sanples,  IJ,  was  chosen  as  10.  This  meant  that 
the  error  in  x  would  be  of  the  order  of  25  per  cent,  vdiich  is 
representative  of  the  errore  in  the  actual  Monte  Carlo  problems  at 
their  maxunun  range  of  usefulness.  If  we  had  chosen  very  large 
sample  sizes,  say  of  the  order  M  ■  1000,  then  the  error  would  be  of 
the  order  of  only  a  few  per  cent.  This  is  a  region  not  particularly 
interesting  from  the  standpoint  of  the  exact  interpretation  of  the 
error.  V/hen  the  error  is,  say,  per  cent  in  a  typical  Monte  Carlo 
problem,  we  are  not  very  often  interested  in  Adiether  it  is  2  per 
cent  or  4  P®r  cent,  and  an  examination  of  the  errors  in  such  regions 
as  this  would  therefore  not  be  particularly  rewarding  for  our 
purposes.  At  the  other  extreme,  if  we  had  taken  samples  of  size  11-2 
or  3,  vdiere  the  error  might  be  in  the  neighborhood  of  75  to  100  per 
cent,  this  too  would  be  a  rather  useless  region  to  explore  since 
the  error  here  would  be  so  large  that  our  knovdng  whether  it  v/as  75 
per  cent  or  125  per  cent  would  be  of  little  use.  A  sample  size  of 
10  was  accordingly  chosen  so  that  the  error  would  be  in  the  range 
where  it  was  most  useful  to  have  a  better  understanding  of  the  exact 


behavior  of  the  error. 
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Having  chosen  II  =  10,  we  then  ran  10,000  runs  of  10  each/’’ 
Histogranis  were  then  obtained  shoiving  the  distribution  functions  for 
X,  o  and  T.  Figures  1,  2  and  3  give  the  results  for  these  three 
distribution  functions  respectivel;^'.  Tables  1,  2  and  3  give  the 
detailed  numerical  results  from  vdiich  these  figures  \^re  prepared. 


•'S'  A  run  of  10  means  that  x,  c 
Eqs.  (6),  (7)  and  (8)  for  II  =  10. 
10,000  tines,  independently. 


and  T  ivere  computed  according  to 
This  procedure  was  then  repeated 
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IV.  RESULTS 

Maxtor  interesting  and  unexjxected  conclusions  can  be  draxvn  by- 
examining  Figs.  1,  2  and  3*  Turning  our  attention  first  of  all  to 
Fig.  1  and  Table  1,  which  give  the  distribution  fvinction  for  x, 
one  notices  that  its  expected  value  is  extremely  close  to  the  actual 
expected  value  of  x,  or  x,  namely  2.0.  However,  the  median  and  most 
probable  values  are  considerably  less  than  the  true  mean  of  2.0.  In 
the  Monte  Carlo  results  obtained  from  the  RAIID  transport  codes,  it 
is  customary  to  draw  by  eye  a  smooth  curve  through  the  various  values 
of  a  histogram.  Even  thoiigh  one  has  drawn  hundreds  of  such  curves, 
it  is  not  ob-vious  which  ones  of  the  above  statistics  are  really 
incorporated  in  such  a  process.  In  other  words,  does  one  try  to 
draw  the  smooth  curve  through  the  expected  value,  the  median  value, 
the  most  probable  value,  or  some  combination  of  these  statistics? 

Mj'  guess  is  that  something  like  the  median  or  the  most  probable 
value  is  more  heavily  weighted  than  the  expected  value.  If  this 
is  true,  then  one  can  see  by  referring  again  to  the  distribution 
shown  in  Fig.  1  that  a  biased  estimate  of  the  mean  is  being  obtadned, 
idiich  in  this  particular  case  can  underestimate  the  true  mean  by 
around  10  per  cent. 

Turning  our  attention  now  to  Fig.  2  and  Table  2,  vdiich  give 
the  distribution  of  the  sample  standard  de-viation  of  the  sample  mean, 
it  can  be  seen  at  a  glance  that  although  this  statistic  has  an 
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expected  value  of  infinity,  a  very  reasonable  sort  of  distribution 
fxmction  is  nevertheless  obtained.  One  can  see,  however,  that  its 
expectancy  tends  to  be  large  because  some  very  large  values  are 
cropping  up  in  the  tail  of  the  curve.  It  is  of  interest  to  note 
that  the  semi-interquartile  range  of  x,  0.?90  is  very  close  to  the 
median  of  o,  0.273,  but  quite  different  from  the  mean  of  o,  0.435. 

A  possible  interpretation  of  these  facts  is  that  the  error  in  x 
will  be  greater  than  o  about  half  of  the  tine  and  less  than  o 
half  of  the  time.  Rather  than  examine  the  distribution  function  for 

o  in  more  detail,  we  will  now  look  at  the  results  for  Student’s 

.  -j:- 

ratio,  which  is  an  error-measuring  statistic  of  more  value  in 
this  case  than  the  standard  deviation  itself.  Referring  now  to 
Fig.  3  iuid  Table  3,  we  see  first  of  all  that  the  distribution 
function  for  T  is  extremely  skew.  The  reason  for  this  lies  in  the 
fact  that  there  is  a  correlation  betvreen  x  and  o.  As  is  well  knowi, 
the  only  distribution  function  where  the  sample  mean  and  the  sample 
variance  are  distributed  independently  is  Gaussitin.  Our  original 
distribution  function  of  Eq.  (3)  is  very  far  from  Gaussian;  hence 
we  expect  a  sizeable  correlation  between  the  sample  mean  and  the 
sa-igDle  variance.  If  the  sample  mean  is  in  error  on  the  low  side, 
then  the  sample  standard  deviation  will  probably  be  a  large 
underestimate.  Student’s  ratio  will  then  be  a  large  negative 
number;  hence  the  e;ctreme  skewness  to  the  left  of  Fig.  3. 

■i:-Student ’  s  ratio  is  a  tern  usually  applied  to  samples  from  a 
normal  parent  population.  However,  here  we  use  the  term  in  the 
more  general  sense,  and  apply  it  to  samples  from  the  distribution 
of  Eq.  (3). 
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In  the  distribution  function  of  Student* s  ratio  for  N  -  oo, 
and  an  original  normal  distribution  one  expects  to  find  about  68 
per  cent  of  the  area  of  the  distribution  curve  between  -1  and  1, 

In  the  Student *s  distribution  for  K  *  10  and  an  original  normal 
distribution  one  finds  some  65  per  cent  of  the  area  between  -1  and  1, 

In  the  distribution  curve  of  Student* s  ratio,  for  our  example,  there 
is  only  48  per  cent  of  the  area  between  -1  and  1.  In  order  to  contain 
68  per  cent  of  the  area  of  Fig.  3  one  must  encorpass  the  range 
-1.8  to  +1.8.  Therefore  in  inteirpreting  the  meaning  of  the  sample 
standard  deviation  in  Monte  Carlo  results  whence  the  probability 
density  function  is  similar  to  Eq.  (3)»  one  should  be  guided  by  the 
nature  of  the  distribution  functions  obtained  in  this  paper  rather 
than  mechanically  using  the  usual  Gaussian  assumption.  A  rather 
complete  discussion  of  the  use  of  Student* s  ratio  for  the  establishment 
of  a  confidence  interval  for  the  mean  is  given  on  pp.  517-519  of  Ref. 


It  seems  likely  that  if  the  assiimed  distribution  function  of 
Eq.  (3)  were  modified  by  increasing  the  ezcponent  of  x  from  3  to  3 
the  same  sort  of  resiiLts  vrould  be  obtained.  In  other  vrords,  even 
tho’ogh  the  variance  would  now  be  finite,  the  distribution  functions 
for  o  and  Student’s  ratio  would  remain  essentially  the  same,  and 
therefore  the  interpretation  of  the  confidence  interval  irauld  be 
unchanged. 
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V.  LIMITATIONS 

One  niust  keep  in  mind  the  assurqstions  which  were  made  in  the 
problem  solved  above.  The  flux  was  assumed  to  be  isotropic  and 
of  one  velocity  (or  energy) .  Also  only  the  case  N  =  10  was  treated. 
Thus  one  cannot  use  the  results  of  this  problem  to  make  quantitative 
predictions  of  confidence  limits  in  actiml  Monte  Carlo  results. 
However  the  qualitative  aspects  of  the  solution  given  here  ivill 
apply  to  Monte  Carlo  results  wiiich  are  obtained  by  the  sampling 
scheme  given  in  the  Appendix. 

Values  of  N  other  than  10  could  easily  be  run  with  the  present 
code  and  more  information  could  thus  be  obtained.  In  view  of  the 
approximate  nat-ure  of  the  assimed  distribution  function,  this  added 
effort  did  not  seem  justified. 
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VI.  COIICLUSIOIIS 


By  means  of  a  sinple  example  we  have  shown  that  the  sample 
variance  in  a  Monte  Carlo  problem  may  be  used  to  estimate  the  error 
of  the  sanple  mean  even  though  the  variance  of  the  distribution 
function  itself  is  infinite,  provided  that  adequate  precautions  are 
taken. 

We  have  also  shown  that  in  the  particular  exajiple  used  a  biased 
estimate  of  the  mean  value  may  be  obtained  v/hich  vriU  lead  to  an 
underestimate  of  the  mean  by  about  10  per  cent. 


13 


■■ 

A  • 

at  S 


o  -  r 
S  S 

—  »<  V 

>  o  S 

fll 

III 


■■■ 

■■■ 

■■■ 

■■■ 

■■■ 

■■■ 

■■■ 

WM 

■HBj 

!■■ 

RIBi 

H 

mmmmmm 

■■■■■■ 

■ 

■ 

■1 

IBBHB 

■■issa 

■■■■■ 

■■■■■ 

E 

ro  -**<>  <3u»nb*jj 


of  sompl*  moan  for  10,000  runt  of  10  toch 


320 


14 


Frequency  distribution  for  sample  standard  deviation  of  sample  mean 
for  10,000  runs  of  10  each 


17 


Table  1 

FREQUEICY  AND  CUMULATIVE  DISTRIBUTIONS  OF  x  FOR  10,000  SAMPLES  OF 

10  EACH 


_  Frequency  Cumilative 

Interval  for  x  Distribution  Distribution 


0 

.1 

0 

0 

.1 

.2 

0 

0 

.2 

.3 

0 

0 

.3 

.4 

0 

0 

.U 

.5 

0 

0 

.5 

.6 

0 

0 

.6 

.7 

0 

0 

.7 

.8 

0 

0 

.8 

.9 

0 

0 

.9 

1.0 

0 

0 

1.0 

1.1 

0 

0 

1.1 

1.2 

26 

26 

1.2 

1.3 

237 

263 

1.3 

1.4 

569 

S32 

1.4 

1.5 

942 

1774 

1.5 

1.6 

1141 

2915 

1.6 

1.7 

1112 

4027 

1.7 

1.8 

1052 

5079 

1.8 

1.9 

865 

5944 

1.9 

2.0 

732 

6676 

2.0 

2.1 

614 

7290 

2.1 

2.2 

488 

7778 

2.2 

2.3 

407 

8IC5 

2.3 

2.4 

277 

8462 

2.4 

2.5 

257 

8719 

2.5 

2.6 

177 

8896 

2.6 

2.7 

148 

9044 

2.7 

2.8 

132 

9176 

2.8 

2.9 

95 

9271 

2.9 

3.0 

89 

9360 

3.0 

3.1 

72 

943^ 

3.1 

3.2 

66 

9498 

3.2 

3.3 

45 

9543 

3.3 

3.4 

53 

9596 

3.4 

3.5 

48 

9644 

3.5 

3.6 

44 

9686 

3.6 

3.7 

23 

9711 

3.7 

3.8 

21 

9732 

3.S 

3.9 

24 

9756 

3.9 

4.0 

24 

9780 

4.0 

4.1 

20 

9800 

_  ^  i=10 

"  =  15  if ,  "i 


Mean  =>  1.99 

Median  =  1.79 

S  eni-int  erquart  ile 
Range  »  0.790 


18 


Table  1 
(Continued) 


Interval 

Frequency 

Cunulat: 

for  X 

Distribution 

Distribi 

4.1 

4.2 

13 

9813 

4*2 

4.3 

17 

9830 

4-3 

4.4 

13 

9843 

4.4 

4.5 

7 

9850 

4.5 

4.6 

16 

9866 

4.6 

4.7 

12 

9878 

4.7 

4.8 

13 

9891 

4.3 

4.9 

5 

9896 

4.9 

5.0 

6 

9902 

5.0 

6.0 

37 

9939 

6.0 

7.0 

26 

9965 

7.0 

8.0 

8 

9973 

8,0 

9.0 

7 

9980 

9.0 

10.0 

5 

9985 

10.0 

11.0 

4 

9989 

11.0 

12.0 

3 

9992 

12.0 

13.0 

1 

9993 

13.0 

14.0 

0 

9993 

14.0 

15.0 

1 

9994 

15.0 

16.0 

2 

9996 

16.0 

17.0 

1 

9997 

17.0 

18.0 

1 

9998 

18.0 

19.0 

1 

9999 

19.0 

20.0 

0 

9999 

20.0 

21.0 

0 

9999 

21.0 

22.0 

0 

9999 

22.0 

23.0 

0 

9999 

23.0 

24.0 

0 

9999 

24.0 

25.0 

0 

9999 

25.0 

100.0 

1 

10000 
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Table  2 

FRBQUEIJCY  AIID  CUIIULATIVE  DISTRIBUTIOIB  OF  o  FOR  10,000  SAIIPLES  OF 

10  EACH 


Inteir/al 

O 

Q  1 

Frequency 

Distribution 

Cumulative 

Distribution 

0 

.01 

0 

0 

1  ^=10  2 

.01 

.02 

0 

0 

°  “Mloo  ^i-  To 

.02 

.03 

3 

3 

V  1=1 

.03 

.04 

8 

n 

.04 

.05 

.05 

.06 

26 

47 

37 

84 

Mean  =  0.435 

.06 

.07 

88 

172 

.07 

.08 

122 

294 

Median  =  .273 

.08 

.09 

159 

453 

.09 

.10 

.11 

.10 

.11 

.12 

162 

2U 

263 

615 

829 

1092 

Seni- int  er  quart  ile 

Ranee  =  .I46 

.12 

.13 

267 

1359 

.13 

.14 

271 

1630 

.14 

.15 

277 

1907 

.15 

.16 

271 

2178 

.16 

.17 

272 

2450 

.17 

.18 

258 

2703 

.18 

.19 

286 

2994 

.19 

.20 

263 

3257 

.20 

.21 

287 

3544 

.21 

.22 

273 

3817 

.22 

.23 

261 

4078 

.23 

.24 

235 

4313 

.24 

.25 

214 

4527 

.25 

.26 

233 

4760 

.26 

.27 

188 

4948 

.27 

.2S 

188 

5136 

.28 

.29 

175 

53n 

.29 

.30 

165 

5476 

.30 

.31 

173 

5649 

•  31 

.32 

168 

5817 

•  3'* 

.33 

160 

5977 

.33 

.34 

169 

6146 

.34 

.35 

130 

6276 

.35 

.36 

143 

6U9 

.36 

.37 

132 

6551 

.37 

.36 

109 

6660 

.38 

.39 

124 

6784 

.39 

.40 

no 

6894 

.40 

.41 

132 

7026 
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Table  2 
(Continued) 


Interval 

for  a 

.41 

.42 

.42 

.43 

.43 

.44 

.44 

.45 

.45 

.46 

.46 

.47 

.47 

.48 

.48 

.49 

.49 

.50 

.50 

.51 

.51 

.52 

.52 

.53 

.53 

.54 

.54 

.55 

.55 

.56 

.56 

.57 

.57 

.58 

.58 

.59 

.59 

.60 

.60 

.61 

.61 

.62 

.62 

.63 

.63 

.64 

.64 

.65 

.65 

,66 

,66 

.67 

.67 

.68 

68 

.69 

69 

.70 

70 

.71 

71 

.72 

72 

.73 

73 

.74 

74 

.75 

75 

.76 

76 

.77 

77 

.78 

78 

.79 

79 

.80 

80 

.81 

81 

.82 

82 

.83 

83 

.84 

84 

.85 

85 

.86 

Frequency- 

Distribution 

104 
80 
87 

105 
67 
90 
62 
67 
74 

55 
61 
61 
72 
66 
47 

56 

57 
54 
54 

42 
38 
37 
40 
37 
47 

43 
25 
32 

34 

35 
29 
31 
21 

19 
^5 
U 

24 
21 
18 

25 

14 

15 
25 
24 

20 


Cionulative 

Distribution 

7130 

7210 

7297 

7402 

7469 

7559 

7621 

7688 

7762 

7820 

7881 

7942 

8014 

8080 

812.7 

6183 

6240 

8294 

8348 

8390 

8428 

8465 

8505 

8542 

8589 

8632 

8657 

8689 

8723 

8758 

8787 

8818 

8839 

8858 

8883 

8697 

8921 

8942 

8960 

8985 

8999 

9014 

9039 

9063 

9083 
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Table  2 
(Continued) 


Frequency 

CuEulati 

Interval 

for  o 

Distribution 

Distribi 

.86 

.67 

24 

9107 

.87 

.88 

16 

9123 

.88 

.89 

15 

9136 

.89 

.90 

19 

9155 

.90 

.91 

22 

9177 

.91 

.92 

19 

9196 

.92 

.93 

13 

92-14 

.93 

.94 

13 

9227 

.94 

.95 

11 

9238 

.95 

.96 

13 

9251 

.96 

.97 

17 

92.68 

.97 

.98 

9 

9277 

.98 

.99 

10 

9287 

.99 

1.00 

9 

9296 

1.00 

1.1 

100 

9396 

1.1 

1.2 

90 

9486 

1.2 

1.3 

69 

9550 

1.3 

1.4 

59 

96U 

1.4 

1.5 

42 

9656 

1.5 

1.6 

31 

9687 

1.6 

1.7 

45 

9732 

1.7 

1.8 

27 

9759 

1.8 

1.9 

22 

9781 

1.9 

2.0 

22 

9803 

2.0 

2.1 

19 

9822 

2.1 

2.2 

11 

9833 

2.2 

2.3 

13 

9846 

2.3 

2.4 

6 

9852 

2.4 

2.5 

11 

9863 

2.5 

2.6 

11 

9874 

2.6 

2.7 

12 

9886 

2.7 

2.8 

3 

9889 

2.8 

2.9 

9 

9698 

2.9 

3.0 

11 

9909 

3.0 

3.1 

6 

9915 

3.1 

3.2 

5 

9920 

3.2 

3.3 

6 

9926 

3.3 

3.4 

1 

9927 

3.4 

3.5 

4 

9931 

3.5 

3.6 

5 

9936 

3.6 

3.7 

3 

9939 

3.7 

3.8 

3 

9942 

3.8 

3.9 

3 

9945 

3.9 

4.0 

1 

9946 
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Table  2 
(Continued) 


Interval 

Frequency 

for  0 

Dist ribut ion 

4.0 

4.1 

3 

4.1 

4.2 

1 

4.2 

4.3 

1 

4.3 

4.4 

1 

4.4 

4.5 

3 

4.5 

4.6 

4 

4.6 

4.7 

1 

4.7 

4.8 

2 

4.8 

4.9 

1 

4.9 

5.0 

4 

5.0 

5.5 

6 

5.5 

6.0 

8 

6.0 

6.5 

2 

6.5 

7.0 

1 

7.0 

7.5 

1 

7.5 

8.0 

0 

8.0 

8.5 

3 

8.5 

9.0 

1 

9.0 

9.5 

1 

9.5 

10.0 

3 

10.0 

15.0 

5 

15.0 

20.0 

1 

20.0 

25.0 

1 

Ciunulative 

Distribution 


9949 

9950 

9951 

9952 
9955 

9959 

9960 

9962 

9963 
9967 
9973 
99S1 

9983 

9984 

9985 
9985 

9988 

9989 

9990 
9993 

9998 

9999 
10000 
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Table  3 

FREQUENCY  AND  CUMULATIVE  DISTRIBUTIONS  OF  STUDENT* S  RATIO,  T,  FOR 

10,000  SAMPLES  OF  10  EACH 


Interval  for  T 


Frequency  Cunulative 

Distribution  Distribution 


-100 

-25 

4 

4 

-25 

-24 

1 

5 

-24 

-23 

0 

5 

-23 

-22 

2 

7 

-22 

-21 

2 

9 

-21 

-20 

2 

11 

-20 

-19 

2 

13 

-19 

-18 

3 

16 

-18 

-17 

6 

22 

-17 

-16 

5 

27 

-16 

-15 

5 

32 

-15 

-14 

14 

46 

-14 

-13 

19 

65 

-13 

-12 

27 

92 

-12 

-11 

29 

121 

-11 

-10 

50 

171 

-10 

-9 

68 

239 

-9 

-8 

101 

340 

-8 

-7 

112 

452 

-7 

-6 

171 

623 

-6 

-5 

316 

939 

-5 

-4.9 

28 

967 

-4.9 

-4.8 

35 

1002 

-4.8 

-4.7 

27 

1029 

-4.7 

-4.6 

40 

1069 

—4.6 

-4.5 

50 

1119 

-4.5 

-4 .4 

34 

1153 

-4.4 

-4.3 

53 

1206 

-4.3 

-4.2 

51 

1257 

-4.2 

-4.1 

43 

1300 

-4.1 

-4.0 

51 

1351 

-4.0 

-3.9 

50 

1401 

-3.9 

-3.8 

53 

1454 

-3.8 

-3.7 

40 

1494 

-3.7 

-3.6 

65 

1559 

-3.6 

-3.5 

57 

1616 

-3.5 

-3.4 

63 

1679 

-3.4 

-3.3 

72 

1751 

-3.3 

-3.2 

72 

1823 

-3.2 

-3.1 

71 

1894 

-3.1 

-3.0 

90 

1984 

-3.0 

-2.9 

104 

2088 

Mean  =  -  1.52 

Median  »  -  0.744 

S  eni-int  erquart il e 
Ran^e  -  1.37 
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Table  3 
(Continued) 


Interval 

for  T 

Frequency- 

Distribution 

Ciunulative 

Distribution 

-2.9 

-2.8 

82 

2170 

-2.S 

-2.7 

106 

2276 

-2.7 

-2.6 

84 

2360 

-2.6 

-2.5 

86 

2446 

-2.5 

-2.4 

101 

2547 

-2.4 

-2.3 

101 

264s 

-2.3 

-2.2 

107 

2755 

-2.2 

-2.1 

111 

2866 

-2.1 

-2.0 

135 

3001 

-2.0 

-1.9 

1^ 

3129 

-1.9 

-1.8 

rL5 

3244 

-1.8 

-1.7 

132 

3376 

-1.7 

-1.6 

143 

3519 

-1.6 

-1.5 

145 

3664 

-1.5 

-1.4 

154 

3818 

-1.4 

-1.3 

156 

3974 

-1.3 

-1.2 

177 

a5i 

-1.2 

-1.1 

171 

4322 

-1.1 

-1.0 

161 

4433 

-1.0 

-0.9 

193 

4676 

-0.9 

-0.8 

205 

4881 

—0.8 

-0.7 

211 

5092 

-0.7 

—0.6 

207 

5299 

-0.6 

-0.5 

219 

5518 

-0.5 

-0.4 

216 

5734 

-0.4 

-0.3 

223 

5957 

-0.3 

-0.2 

224 

6181 

-0.2 

-0.1 

240 

6421 

-0.1 

0 

255 

6676 

0 

0.1 

266 

6942 

0.1 

0.2 

26? 

7204 

0.2 

0.3 

273 

7482 

0.3 

0.4 

275 

7757 

0.4 

0.5 

294 

8051 

0.5 

0.6 

267 

83I8 

0.6 

0.7 

277 

8595 

0.'' 

0.8 

255 

8850 

0.8 

0.9 

268 

9118 

0.9 

1.0 

^40 

9358 

1.0 

1.1 

192 

9550 

1.1 

1.2 

136 

9686 

1.2 

1.3 

94 

9780 

1.3 

1.4 

72 

9852 

Table  3 
(Continued) 


Frequency  Cumulative 

Interval  for  T  Distribution  Distribution 


1.4 

1.5 

59 

99n 

1.5 

1.6 

38 

9949 

1.6 

1.7 

15 

9964 

1.7 

1.8 

13 

9977 

1.8 

1.9 

10 

9987 

1.9 

2.0 

5 

9992 

2.0 

2.1 

1 

9993 

2.1 

2.2 

3 

9996 

2.2 

2.3 

1 

9997 

2.3 

2.4 

1 

9998 

2.4 

2.5 

0 

9998 

2.5 

2.6 

1 

9999 

2.6 

2.7 

1 

10000 
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Appendix 

DERIVATIOII  OF  BASIC  PROBABILITY  DEMSITY  FUNCTION  FOR  TIE  FLUX 

There  are  several  ways  of  recording  flux  or  some  function  of 
the  flux,  such  as  dosage,  in  a  Monte  Carlo  code.  One  common  scheme 
used  in  the  RAND  neutron  and  gamma  ray  transport  codes,  as  well  as 

(■3) 

in  codes  developed  by  Los  Alamos  and  Sandia  Corporation,  is  to 
record  certain  functions  of  a  particle  as  it  passes  through  a  given 
plane  in  the  sar^ile  space.  The  parameters  which  need  be  kno\i/n  are 
particularly  the  statistical  \veight  of  the  particle,  its  energy,  and 
the  angle  of  crossing  with  respect  to  a  normal  to  the  plane. 

Generally  speaking,  particles  %7ill  cross  any  one  point  on  the 
plane  with  many  values  of  statistical  weight,  energy,  and  angle. 

Here  we  ;vill  consider  the  simplest  case  vdiich  appro;cimates  these 
corplex  conditions.  In  particular,  vje  assume  that  all  particles 
have  the  same  weight,  the  same  energy,  and  that  their  directions 
are  distributed  isotropicsilly  in  space. 

Let  0  be  the  angle  of  a  particle  with  respect  to  the  normal 
to  the  recording  plane.  The  probability  density  function  for  a 
particle  at  angle  0  is  simply 

dp  =  I  sin9  de  0  <  e  <  u.  (Al) 

If  we  consider  particles  vdiich  cross  a  fixed  area  of  the  recording 
plane,  then  the  probability  density  function  for  a  particle  crossing 
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this  area  of  the  plane  at  angle  Q  is 


dp  «  sine  |cose|  de  0  <  0  <  ir. 


(A2) 


In  the  Monte  Carlo  code,  vhen  a  particle  actually  crosses  the 
recording  plane  at  an  angle  9,  the  flux  or  dosage  contributed  by  that 
particle  is  then  given 


cp 


(A3) 


aside  from  a  constant  of  proportionality,  which  we  may  neglect.  If 
vre  eliminate  0  from  Eqs.  (2)  and  (3),  we  obtain  the  probability 
density  function  for  the  flux  per  particle  crossing,  which  is 

ap  ”  9  =  1  to  00  (A4) 


“0  9  “  0  to  1. 

It  is  this  equation  which  the  \ifhole  substance  of  this  mencrandum  is 
based  upon,  Ifcte  that  the  expected  value  of  the  flux  in  Eq,  (4) 
is  2,  This  sirply  means  that  the  average  flux  of  particles  crossing 
a  plane  in  an  isotropic  field  is  twice  the  value  vdiich  would  be 
obtained  by  having  the  same  number  of  particles  cross  at  normal 
incidence,  and  this  is  in  fact  a  \irell-known  result. 
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ADDENDUM  ON  RUMIMG  TIME 


The  running  time,  including  printing,  for  the  10,000  san?iLes 
of  10  each  was  4^  minutes  on  the  7090, 
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